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Abstract 

The classical EMD algorithm has been used extensively in the literature to decom- 
pose signals that contain nonlinear waves. However when a signal contain two or more 
frequencies that are close to one another the decomposition might fail. In this paper 
we propose a new formulation of this algorithm which is based on the zero crossings of 
the signal and show that it performs well even when the classical algorithm fail. We 
address also the filtering properties and convergence rate of the new algorithm versus 
the classical EMD algorithm. 
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1 Introduction 



In scientific literature there exist many classical sets of functions which can decompose a 
signal in terms of "simple" functions. For example Taylor or Fourier expansions are used 
routinely in scientific and engineering applications (and many other exist). However in all 
these expansions the underlying functions are not intrinsic to the signal itself and a precise 
approximation to the original signal might require a large number of terms. This problem 
become even more acute when the signal is non-stationary and the process it represents is 
nonlinear. 

To overcome this problem many researchers used in the past the "principal component 
algorithm" (PCA) to come up with an "adaptive" set of functions which approximate a 
given signal. A new approach to this problem emerged in the late 1990's when a NASA 
team has developed the "Empirical Mode Decomposition" algorithm(EMD) which attempts 
to decompose a signal in terms of it "intrinsic mode functions" (IMF) through "sifting 
algorithm". A patent for this algorithm has been issued [1]. 

The EMD algorithm is based on the following quote [2]: "According to Drazin the first 
step of data analysis is to examine the data by eye. From this examination, one can immedi- 
ately identify the different scales directly in two ways: by the time lapse between successive 
alterations of local maxima and minima and by the time lapse between the successive zero 
crossings.... We have decided to to adopt the time lapse between successive extrema as the 
definition of the time scale for the intrinsic oscillatory mode" 

A step by step description of the EMD sifting algorithm is as follows: 

1. Let be given a function f(t) which is sampled at discrete times k = 1, . . . n}. 

2. let h (k) = f(t k ). 

3. Identify the max and min of ho(k). 

4. Create the cubic spline curve M x that connects the maxima points. Do the same for 
the minima M n . This creates an envelope for h (k). 

5. At each time t k evaluate the mean m k of M x and M n {m k is referred to as the sifting 
function). 

6. Evaluate h±(k) = h (k) — m k . 
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7. If the norm of \\ho — hi\\ < e for some predetermined e set the first intrinsic mode 
function IMF X = hi (and stop). 

8. If the criteria of (7) are not satisfied set h (k) = hi(k) and return to (3) ("Sifting 
process"). 

The algorithm has been applied successfully in various physical applications [1-6]. How- 
ever as has been observed by Flandrin [3] and others the EMD algorithm fails in many cases 
where the data contains two or more frequencies which are close to each other. 

To overcome this difficulty we propose hereby a modification of the EMD algorithm by 
replacing steps 4 and 5 in the description above by the following: 

4. find the midpoints between two consecutive maxima and minima and let Nk be the 
values of ho at these points. 

5. Create the spline curve that connects the points N^. 

The essence of this modification is the replacement of the mean which is evaluated by 
the EMD algorithm as the average of the max-min envelopes by the spline curve of the 
mid-points between the maxima and minima. This is in line with the observation by Drazin 
(which was referred to above) that the scales inherent to the data can be educed either 
from the max-min or its zero crossing. In the algorithm we propose hereby we mimic the 
"zero-crossings" by the mid-points between the max-min. 

It is our objective in this paper to justify this modification of the EMD algorithm through 
some theoretical work and case studies. The plan of the paper is as follows: In Sec 2 we 
provide theoretical justification for the new algorithm proving that it acts as a high pass 
filter for certain classes of signals. In Sec. 3 we provide examples of a signal composed of 
two or three close frequencies (with and without noise) where the classical EMD algorithm 
fails but the modified one yields satisfactory result. In Sec. 4 we discuss the convergence 
rate, resolution and related issues concerning the classical and new "midpoint algorithm". 
We end up with some conclusions in Sec 5. 

2 Theoretical Justification 

In this section we provide a theoretical justification for the proposed modified EMD algorithm 
by analyzing it performance on several generic signals which contain several close frequencies. 
However in this analysis linear, quadratic and cubic interpolating polynomials will be applied 
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to represent the midpoints interpolating function (instead of splines). To motivate this 
"replacement" we observe that the coefficients of each spline polynomial depend non- locally 
on the data i.e. these coefficients might change if additional data is added. On the other 
hand Lagrange interpolating polynomials depend only the local data. 
Lemma 1: Consider a signal of the form 

f(t) = f 1 (t) + f 2 (t) + f 3 (t) (2.1) 

where 

ft(t) = cos(ut), f 2 {t) = cos((l + ae)ut), f 3 (t) = cos((l + be)ut), b>a>0, < e < 1. 

(2.2) 

Let the projection of the midpoint linear interpolating function on fi(t), fift), fs(t) over 
an interval containing five midpoints be denoted respectively by Pu, P 2 i, P31 then 

3 UJ 

and 

Pu _ P3i ^(l±^j^ e3+0(e% (23) 

O UJ 

P„- fil » '^°'- h )V + (^ (2.4) 

3 UJ 

p il -P 31 = 2 ^» 2 + a2 -'"'^- a K^O^). (2.5) 
3 uj 

Proof: As a first step we find the approximate location of the extrema of f(t) on the 
interval [0,67r]. To do so we differentiate f(t) and observe that due to the fact that e C 1 
the locations of these points are close to — . Setting t = — + rj we expand f'(t) in a Taylor 
series to order 2 in 77 around — . Then we solve for 77 to obtain the approximate locations 
of the extrema points. Taking the five midpoints {t\, . . . ,t$} between these extrema and 
evaluating f(ti), i = 1, ... ,5 we construct the linear interpolating function gi(t) between 
these points. The projection of fi(t) on g\{t) is 

Pu = / fi{t)gi{t)dt. (2.6) 



Expanding Pu — Pji in a Taylor series in e one obtains (I2.3p - fl2.5l) . 
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Lemma 2: With the same settings as in Lemma 1 let the function Q2(t) consists of the two 
quadratic polynomials interpolating {ti,t 2 ,t 3 } and {£3, £4, £5} respectively. The differences 
between the projections 

Pi2= I " fi(t)g 2 (t)dt (2.7) 
Jti 

are given by eqs (12.3p . (12.4p and (12.51) respectively (where Pn is replaced by P^) 

Proof: As in Lemma 1 we compute the projections Pa, and expand the results in a 

Taylor series in e to obtain (I2.3l) - (l2.5p 

Lemma 3: With the same settings as in Lemma 1 if the projection of fi{t), f2{t), f3{t) is 

made on the cubic polynomial g 3 (t) interpolating {£1, £2, £3, £4} then the differences between 

the projections 

Pa= I* fi(t)g 3 (t)dt (2.8) 

are given by 



ti 



Pl3 -P S ^ + ° 2 -^-^ W). (2.9) 
9 con 

P 13 _P 23 = !^+° 2 -M(49^-57) e 3 + 0(ea); (210) 
9 con 

9 LU7V 

Theorem 1: As a result of one iteration of the midpoint EMD algorithm with linear, 
quadratic or cubic interpolating functions the change in the projections of the functions 
fi(t) i = 1, 2, 3 on the signal in the interval [ti, £5]) ([£1, £4] in the cubic case) satisfy AAij > 
AA 2) j > AA 3 j, j = 1,2,3. (Here j represents the different interpolating functions). 

Proof: The projection of fi(t) on the original signal is 

4, = f f(t)fi(t)dt. 
After one iteration the signal is represented by 

f\t) = f(t) - 9j (t) 
and the projection of fi(t) on f 1 ^) is 

4, = /*(/(*) - 9j(t))fi(t)dt = A\ 3 - P id . 
Ju 
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Hence 

A-l, = 4i - A\ 3 /',, >0 

From the results of lemmas 1, 2, 3 we have that For j = 1, 2, 3 Py — P 2 j > 0, Py — P37 > 
and P 2i - P 3j > 0. It follows then that AAy > AA 2 j > AAjj > 0. 

We conclude therefore that in the new signal (after one iteration) the amplitude of fs(t) 
will be larger than those of f2{t) and fi(t). In other words the midpoint EMD algorithm 
acts as a high pass filter. 

We consider now a signal with two close frequencies where a phase shifts exists between 
these two frequencies. 

Lemma 4: Consider a signal of the form 

f(t) = h(t) + f B (t) (2.12) 

where 

f 4 (t) = cos(wt), f s (t) = cos[(l + ae)wt + <f>]. (2.13) 

where a > 0, and < e,^ C 1. With same setting as in Lemmal let the projection of 
the midpoint linear interpolating function (for f(t) defined in (j2.13|) ) on /^(t), fs(t) over an 
interval containing five midpoints {£i, . . . ,t$} be denoted respectively by P41, P51 then 

_ = Mfaw + ts^ + ^n + 4 _ ^ 4 

Lemma 5:With the same settings as in LemmaA let the function g^{t) consists of the two 
quadratic polynomials interpolating {^1,^2,^3} and {£3, £4, £5} respectively. The difference 
between the projections P 42 , P52 of f^(t), f§{t) on g 5 (t) is given by f)2.14p (where Pn is 
replaced by P i2 ). 

Lemma 6: With the same settings as in Lemma 4 let the function ge(t) consists of the 
cubic polynomial interpolating {ti, t 2 , t$, £4}. The difference between the projections P43, P53 
of U(t), f 5 (t) on g 6 (t) is 

p 43 _p 53 = J^— r a 2 vr 2 (497r 2 - 57)e 3 + 8a7T(/)(57r 2 - 6)e 2 + 20 2 (5^ 2 - 6)el +0(e 4 , 4 ) (2.15) 

Theorem 2: As a result of one iteration of the midpoint EMD algorithm with linear, 
quadratic or cubic interpolating functions the amplitudes B it j, i = 1,2 j = 1,2, 3 (where the 
index j represents the different interpolating functions) of the two frequencies present in the 
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signal (I2.12p - |2~T3|) ) will satisfy Bij < B2J. In other words the midpoint EMD algorithm for 
this signal is a high pass filter. 

Proof: The proof is similar to the proof of theorem 1. 

2.1 Perturbation Analysis 

To investigate the performance of the EMD algorithm (classical and midpoint) in the presence 
of a perturbation (viz. noise) we considered a signal of the form 

S (t) = cos(ut) + ef(t) (2.16) 

where 0<e<l. To analyze this signal we assume that the presence of noise (represented 
by ef(t)) does not change (appreciably) the location of the extrema in the signal i.e. the 
maximum and minimum are located respectively at the following times 

2/C7T (2k+l)lT , 

Pk = , Qk = - fc = 0,l,.... 2.17 

UJ UJ 

The value of the signal at these points is 

S (p h ) = l + ef(—), S (q k ) = -l + ef( ( - 2k+1)7r ), k = 0,l,.... (2.18) 

UJ UJ 

To apply the classical EMD algorithm to this data one has to compute the spline curves 
Smaxif) and S min (t) for the points (pk,S(Pk)) an d (qk,S(<lk)) respectively. The new signal 
after one iteration of the (classical) EMD algorithm is given by 

S°{t) = S (t) - Smax ^ I Smm ^ (2.19) 

Similarly for the new EMD algorithm we take the midpoints dj = ^^7^-, j = 0,1,... 
between the extrema of the signal and evaluate the signal at these points to obtain 

So(d j ) = ef(W±pl). (2.20) 

Computing the spline curve S mi d for the data points (dj, S (dj)), and subtracting this from 
the original signal we obtain after one iteration of this algorithm that the new signal is given 
by 

S?(t) = S (t)-S md (t). (2.21) 
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To compare the noise reduction efficiency of the two algorithms for this signal on a finite 
time interval [0, ( 2n+1 ) 7r ] (i. e k — 0, . . . ,n and j = 0, . . . , 2n) we project the new signals on 
{cos(a;t), sin(u;t)}. (Both {cos(ut) , sm(ut)} have to be considered due to a possible phase 
shift in the new signal). To this end we have to compute 

P{= \ S\{t) cos(cut)dt, P 2 C = / S$(t) sm(oot)dt. (2.22) 

Jqo J go 

Ql= I S?(t)coa(ut)dt, Q n 2 = \ S?(t)sm(ujt)dt. (2.23) 

Jdo Jdo 

Using fl2TT6l) - fl2~2T|) yields 

PI = ^ + J Pn [ef{t) - S ™*W + S ™"M ] cos ( ut )dt (2.24) 

P - = [ € f(t) - S max( t ) + S min(t) ] ^.25) 

Jqo 2 



d'2r, 



Qi = 7T+ f e /(*) - S mid(t)] cos(cot)dt (2.26) 

ZUJ Jdo 

Q2 = / Wit) - S mU {t)\ sm(oot)dt. (2.27) 

Jdo 

We conclude then that the efficiency of the algorithm to eliminate the noise in the signal can 
be measured by the smallness of the absolute values of the integrals 

fc [f/(() _U) + Ut) ]cos(u() ^ Qmid= j d2 \ tf{t) _s rmd {t)]^{ujt)dt 

<?0 ^ Jdo 

and the absolute values of P£, Q%- 

To obtain a quantitative insight into this issue we considered the special case where 

f(t) = cos(^). 

with k = 0, . . . , 9 and j = 0, . . . , 18. A calculation of the P mn and the other integrals for 
v 00 yields: 

26.703e 13.3526, . x2x 

Pmn = + 5— (f - W) + 0((U - uf) 

OJ or 

28.274e 42.41e. . , 2 . 

Qmid = + o— -U)) + 0((U - Uf) 

OJ 0J A 

p 2 «==^(,-.) +0 ((,-^) 



Ql = '—{v -u>) + 0((u - uf) 



-12.207e 
~ 2 

These results show that when the frequency of the noise is close to the original frequency 
the classical algorithm leads to a large phase shift in the signal and the noise is shifted with 
it. 

In a more general setting of this analysis one may consider a Fourier expansion of f(t) if 
this function is periodic. 

For the convergence of the sifting iteration we now prove the following: 

Theorem 3: For the signal ( 12 . 1 6[) if we replace the spline approximation between the 
midpoints by a linear interpolating function then the sifting process will converge to cos(uit) 
if the derivatives of f(t) in the L 1 norm are bounded. 

Proof: The coordinates of the the midpoints between the max-min of (12.161) are 

d 3 = (—2^—, = 1, + 1 



For a linear interpolation function eq. (12 .2 1 j) becomes 



. x , u [A— 2£— J -JK—to-) V —) ^(2^ + 1^, 

S?(t) = cos(ut) + e { f{t) - ^ J /( V ^ ^ ; 

j=o 7r 

(2.28) 

The L 1 norm of the signal So(t) is 0(1). To obtain an approximation for the L 1 norm 
of the perturbation Pi = ef(t) — S^ id after one iteration we use trapezoidal integration. In 
this setting the integral of f(t) cancels the integral of the linear interpolating function. This 
yields the following standard estimate for the residue of the perturbation 

\\Pi\\ = eO((-r)\\f"(t)\\ 

(where primes denote differentiation with respect to t). We conclude that if the L l norm of 
derivatives of f(t) are bounded then the sifting iterations will converge. 

Using the same settings as in theorem?) (i.e replacing the spline interpolating function 
by a linear interpolating function) similar results apply to an d Q2 

Lemma 7: In the L l norm 

110? " ?ll = zO((*f)\\f"(t)\l \\CK\\ = eO((-n\f"(t)\\ 
The proof is the same as in theorem 3. 



(2j + l)7T- 
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2.2 Some Additional Analytical Insights 

To obtain analytical insights about the performance of the EMD-midpoint algorithm we 
considered a signal of the form 

f(t) = ^[cos(wit) + cos(w2*)], (2.29) 

where the ratio of the frequencies ui2 is a rational number viz. 

u 2 m 
u>i n 

where m,n are relative prime integers. In this case the signal f(t) is actually periodic with 
period p = Due to this fact behavior of the classical versus the mid-point algorithm can 
be delineated without the need to discretize the signal. 

On the interval [0,p] the extrema of the signal which satisfy ^ = are given by 

sin Uit ui2 m 
sin u 2 t ui\ n 

Computing these extrema points it is straightforward to construct the spline approximations 
Smax(t), S min {t) to the maximum and minimum points and compute their average. Similarly 
we can find the midpoints between the maxima and minima and evaluate the corresponding 
spline approximation S m id{i) to the signal at these points. After one iteration of the sifting 
process the "sifted signal" is given respectively by 

U) = /(*) - Smax{t) + Smm{t \ (2.30) 

and 

hmid(t) = f(t) - S mid (t). (2.31) 

The efficiency of the two algorithms can be deduced by projecting these new signals on the 

Fourier components of the original signal. To this end we compute 

rp rp 
a mn = / h mn (t)cos(u4)dt, b mn = / h mn (t) s\n(u A t)dt. (2.32) 
Jo Jo 

rp rp 

Cmn= / h mn (t)cos(uj 5 t)dt, d mn = / h mn (t)sm(u) 5 t)dt. (2.33) 

Jo Jo 

and 

rp rp 

a m id= / h mid (t)cos(uj4t)dt, b m i d = / h mid (t) sm(u} 5 t)dt. (2.34) 

Jo Jo 
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rv rp 
Cmid = / h mid (t) cos(u 4 t)dt, d mid = / h mid (t) sxn(u 5 t)dt. (2.35) 
Jo Jo 

The amplitude of the Fourier components of the two frequencies in the classical EMD algo- 
rithm is 

A mn = \/ a ran ~i~ ^mn' Bmn = \/ c mn ^mri' (2.36) 

Similarly for the mid-point algorithm we 

Amid = ytfnid + Kiid> B rnid = \J tfnid + ^mid' (2.37) 

The objective of the sifting process is to eliminate one of the Fourier components in favor 
of the other. As a result the first IMF will contains, upon convergence, only one of the 
Fourier components in the original signal. Therefore the efficiency of the two algorithm can 
be inferred by comparing A mn versus B mn and A mid versus B mid . 
In the particular case where the signal is given by 

1 3tT 7T 

/(*) = ^[cosi^t) + cos(w2*)], wi = -rr, ^2 = — • (2.38) 
p = 128 (See Fig. 5). Computing the integrals that appear in eqs. f l2.32p - fl2.35p we obtain 

A™ = 31.63346911, B mn = 29.70292046, (2.39) 

A mid = 34.19647843, B mid = 20.81145369. (2.40) 

These results show that after one iteration the classical EMD did not separate the two 
frequencies effectively (A mn and B mn are close to each other). On the other hand the mid- 
point algorithm performed well. 



3 Examples and Comparisons 

Extensive numerical experiments were made to test and verify the efficiency of the modified 
algorithm. We present here the results of one of these tests in which the signal contains 
three close frequencies (where the classical EMD algorithm fails). In our tests we considered 
also the effects of noise and phase shifts among the different frequencies but these will not 
presented here. 

/(£) = -[cos(wit) + cos(u2t) + cos(a>3t)] (3.1) 
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where 

7T 

wi = 12cj , uj 2 = lOwo, w 3 = 8w , = 77777- 

25o 

To apply the new EMD algorithm to this signal, discretized it over the time interval [—2048, 2048] 
by letting t k+1 - t k = 1, k = 1, . . . , 4097. 

The results of the signal decompositions into IMFs are presented in figures 1—4. In 
all these figures the red lines represent the frequencies in the original signal (or its power 
spectrum) and the blue lines the corresponding intrinsic mode functions or their power 
spectrum which were obtained by the midpoint algorithm. 

Fig. 1 is a plot of the data for the signal described by (13. ip . Fig. 2 represents the first 
IMF in the decomposition (versus the leading frequency in the data) while Figs. 3 — 4 depict 
the spectral density distribution for the first two IMFs versus those related to the original 
frequencies in the data. It should be observed that although the amplitude of the spectral 
densities in these plots are somewhat different the maxima of the spectral density in each 
plot is very close to the original one. 



3.1 Cubic Lagrange Interpolation 

In both classical and the new versions of the EMD algorithm splines are used for interpolation 
purposes. However the coefficients of each spline polynomial depend non-locally on the data. 
As a result these coefficients might change if additional data is added. To compare the two 
algorithm without this non-local dependence we replaced the spline interpolation by cubic 
Lagrange interpolation (where the coefficients of the interpolating polynomial depend on the 
local values of the function on the interval). 

To carry out this comparison between the two EMD algorithms we considered a signal 
composed of two frequencies and noise, 

3 

S (t) = cos(ujt) +cos(-cot) + ef(t), e < 1, (3.2) 

on the time interval [0, ^]. The time interval was chosen so that the signal (without the 
perturbation) has four maxima and minima on this interval. 

As in subsection 2.1, we assume that the locations of the maxima and minima do not 
change appreciably due to the perturbation. These locations are then given respectively by 



2 / V25-2VT0 \\ 4tt 4vr 

p = 0, pi = - 7r - arctan , p 2 = p%, p3 = — , (3.3) 

u \ \ 1 + a/10 / / oj u) 
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q = - arctan =- I gi = — , g 2 = g , g 3 = — + q . (3.4) 

w y \ _1 + vlO J J co u co 

Computing the cubic Lagrange interpolating polynomials L max and L min for the maxima 
and minima respectively we obtain for the "modified Lagrange classical-EMD" after one 
iteration 

Sfc) = S (t) - Lmaxit) + Lmm(t) . (3.5) 

The number of midpoints on the interval [0, ^] is seven. For this reason we use two cubic 
Lagrange interpolating polynomials on this interval. (The first is valid over the interval 
[c?o,c?3] and the second is valid over [c^,^]). Denoting this combined polynomial by L mic i 
and subtracting from the original signal we obtain after one iteration of this algorithm that 
the new signal is given by 

S?(t) = S (t)-L mid (t). (3.6) 

To examine the performance of the two algorithms we project these new signals on cos(|w) 
and sin(|u;). 

pc= / S1(t) cos(-Lot)dt, P 2 C = / S1(t)sm(-ujt)dt. (3.7) 

J q 2 J qo 2 

aIld A ,1 

Q{= \ S?(t)cos(-ut)dt, Q n 2 = \ S?(t)sm(-iot)dt. (3.8) 
Jdo 2 J do 2 

Furthermore if we assume that / = cos(z4) and v ~ \oo we obtain to order e 

^ 3.8568 0.0175e ^. 3 , nr 1.0637 0.0399e 3 , 

Pr = + +0(i/--w), P 2 6 = + + 0(u--u) (3.9) 

WW Z CO co z 

^ n 6.3795 0.1257e 3 , ^ n 0.2184 0.3113e 3 , . oin , 

Qi = + O(v--co), Ql = + + O(u--oj). (3.10) 

CO LO Z CO CO Z 

These results demonstrate the superiority on the midpoint algorithm in this setting. (The 
total projection of the new signal on cos §o>t is larger and the phase shift is smaller). 



4 Convergence Rates 

To compare the convergence rates of the classical versus the midpoint algorithm we consid- 
ered three cases all of which were composed of two frequencies. In the first case the two 
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frequencies were well separated. In the second case the two frequencies were close while in 
the third case they were almost "overlapping". In all cases the signal was given by 

f(t) = i(cosojit + coso> 2 £) 

This signal was discretized on the time interval [—2048,2048] with At = 1. 
For the first case the two frequencies were 

71 

cji = 12a;, o>2 = 8a;, u 
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As can be expected both the classical and midpoint algorithm were able to discern the indi- 
vidual frequencies through the sifting algorithm. However it took the classical algorithm 59 
iterations to converge to the first IMF. On the other hand the midpoint algorithm converged 
in only 7 iterations (using the same convergence criteria). We wish to point out also that the 
midpoint algorithm has a lower computational cost than the classical algorithm. It requires 
in each iteration the computation of only one spline interpolating polynomial. On the other 
hand the classical algorithm requires two such polynomials, one for the maximum points and 
one for the minimum points. 

For the second test the frequencies were 

71 71 7Y 71 
Wi = 1 , U)o = 

1 24 288' 24 288 

that is the difference between the two frequencies is j^. 

In this case the midpoint algorithm was able to separate the two frequencies. Fig 6 and 
Fig 7 compare the power spectrum of the original frequencies versus those of IMFi and 
IMF2 which were obtained through this algorithm. Convergence to IMFi was obtained in 
18 iterations and IMF 2 was obtained by 7 additional iterations. 

The classical EMD algorithm did converge to IMFi in 45 iterations but the power spec- 
trum of this IMF deviated significantly from the first frequency in the signal. IMF 2 failed 
(completely) to detect correctly the second frequency. 

In third case the frequencies were 

u l — 77T + 7777777) u l 



24 1000 24 1000 

In this case the classical algorithm was unable to separate the two frequencies i.e IMFi 
contained both frequencies. The midpoint algorithm did somewhat better but the resolution 
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was not complete. Moreover the sifting process in both cases led to the creation of "ghost 
frequencies" which were not present in the original signal. 

At this juncture one might wonder if a "hybrid algorithm" whereby the sifting function 
is the average (or some similar combination) of those obtained by the classical and midpoint 
algorithms might outperform the separate algorithms (in spite of the obvious additional 
computational cost). However our experimentations with such algorithm did not yield the 
desired results (i.e. the convergence rate and resolution did not improve). 

5 Conclusion 

In this paper we presented a variant of the EMD algorithm which utilizes the midpoints 
between the max-min points of the signal in the sifting iterative process. We demonstrated 
through several case studies and theoretical approximations that this algorithm can resolve 
signals with moderately close frequencies where the classical EMD algorithm fails. We 
showed also that it has a better convergence rate. From a formal point of view this superior 
performance of the midpoint algorithm can be traced to the fact that the deviation of the 
signal average from zero is sampled at "half the scale of the classical EMD algorithm. 

References 

1 N. E. Huang - USA Patent #6,311, 13051 , Date Oct 30,2001 

2 N. E. Huang et all, "The empirical mode decomposition and the Hilbert spectrum for 
nonlinear and non- stationary time series analysis", Proceedings of the Royal Society 
Vol. 454 pp. 903-995 (1998) 

3 Gabriel Rilling and Patrick Flandrin, "One or Two Frequencies? The Empirical Mode 
Decomposition Answers", IEEE Trans. Signal Analysis Vol. 56 pp. 85-95 (2008). 

4 Zhaohua Wu and Norden E. Huang, "On the Filtering Properties of the Empirical 
Mode Decomposition, Advances in Adaptive Data Analysis", Volume: 2, Issue: 4 pp. 
397-414. (2010) 



15 



5 Albert Ayenu-Prah and Nii Attoh-Okine, "A Criterion for Selecting Relevant Intrin- 
sic Mode Functions in Empirical Mode Decomposition", Advances in Adaptive Data 
Analysis, Vol. 2, Issue: 1(2010) pp. 1-24. 

6 G. Rilling, P. Flandrin and P. Goncalves, "Empirical Mode Decomposition As a Filter 
Bank, IEEE Signal Processing Letters, vol. 11, no. 2, pp. 112-114, 2004 



16 



IMF 1 vs. first frequency in the data 
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Figure 5: 
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Power Spectrum of IMF 1 vs. first frequency 
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Power Spectrum of IMF 2 vs. 2nd frequency 
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